The dynamics of condensate shells: collective modes and expansion 
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We explore the physics of three-dimensional shell-shaped condensates, relevant to cold atoms in 
"bubble traps" and to Mott insulator-superfluid systems in optical lattices. We study the ground 
state of the condensate wavefunction, spherically-symmetric collective modes, and expansion prop- 
erties of such a shell using a combination of analytical and numerical techniques. We find two 
breathing-type modes with frequencies that are distinct from that of the filled spherical conden- 
sate. Upon trap release and subsequent expansion, we find that the system displays self-interference 
fringes. We estimate characteristic time scales, degree of mass accumulation, three-body loss, and 
kinetic energy release during expansion for a typical system of *^Rb. 
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The cooling and trapping of dilute atoms has recently 
achieved unprecedented levels of control and sophistica- 
tion. With the advent of optical lattices 111, quasi-one 
and -two dimensional trapping potentials [3, and mix- 
tures of different species [3| , condensates of bosonic atoms 
have been created in a plethora of interesting geometries. 
Boson mixtures in a particular regime of interactions can 
form a layered or core-and-shell structure |j]. Bosons 
in a three-dimensional optical lattice also display a shell 
structure as a result of the confining trap |3- It may 
even be possible to confine a dilute atomic condensate to 
a spherical shell-shaped region by means of a specifically- 
designed "bubble trap" |g. Towards an understanding 
of these, and other systems where superfluid order ex- 
ists in "hollow" geometries, we consider the physics of a 
condensate whose shape is a three-dimensional spherical 
shell. We identify key features in the condensate collec- 
tive modes and expansion upon trap release that would 
distinguish such shell-shaped condensates from the more 
common filled cases. Moreover, we find that expansion 
properties have distinct similarities with the "Bose-nova" 
experiments \j\. 

As a potentially striking application, our study of con- 
densate dynamics ought to provide a novel probe for 
the predicted co-existence of shell-shaped superfluid and 
Mott insulator regions in optical lattices [S| ■ At the mean 
field level, and at scales much larger than the optical lat- 
tice spacing, the physics of the superfluid regions should 
be well-described by a smooth effective potential repre- 
senting the Mott regions. In fact, our analyses ought 
to be relevant to any mean-field description which yields 
shell-shaped condensate confinement, as for instance in 
the case of boson mixtures where repulsive interactions of 
one species exclude another species from occupying the 
center of the trap. 

We therefore consider the ground state wavefunc- 
tion for the condensate assuming an effective three 
dimensional confining potential of the form V(y) = 
(l/2)TOijjQ(r — ro)^ (where m is the atomic mass). The 
ground state wavefunction of the condensate obeys the 



time-independent Gross-Pitaevskii (GP) equation: 
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where 5 is a measure of the repulsive interactions between 
the atoms ((7 = ATrfi^a/m with a the s-wave scattering 
length) and /x is the chemical potential, set by the nor- 
malization condition: J |'0(r)| (Pr = N. The density of 
condensed atoms is given by n(r) = |V'(r)p. 

For suitably large N, the kinetic energy is very small 
compared to the potential energies and the gradient 
term may be neglected (the Thomas-Fermi approxima- 
tion) [SJ, giving an approximate ground state wave- 
function xPtf{y) = yj[^ — V{r)\/g in the regions where 
/i > V{v) and zero elsewhere. In this approximation, 
the condensate occupies a spherically-symmetric shell 
centered at rp with half-width ri = y^2fi/{mu!Q). In 
the limit of a thin shell (ri <^ ro), ri is found to be 
[3gN / {SnmujQr'^)]^^^ . Therefore, for the Thomas-Fermi 
wavefunction, the condition for a thin shell can be writ- 
ten 
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where we have introduced the oscillator length Oosc = 
y'h/JmujQ). In order to obtain analytic results, we will 
often work within this "thin shell" limit. While appli- 
cable to typical conditions in bubble traps and optical 
lattice systems, we also expect an analysis of thin shells 
to capture the salient features of thicker shells. Notable 
expected differences will be discussed in some cases. 

While we are able to go beyond the Thomas-Fermi 
approximation using numerical techniques, it is worth 
estimating its regime of validity. Using a radially sym- 
metric Gaussian wavefunction centered at radius Rg and 
with characteristic width i?i, we find the ratio of the 
zero-point kinetic energy to the potential energy of in- 
teraction to scale as i?o/(iVai?i). (The analogous ratio 
for a condensate in a harmonic trap centered at r = 
scales as R/{Na) (see, e.g., |j|)). Hence, for the shell 



(i.e., taking Rq to be ro and i?i to be ri), the Thomas- 
Fermi approximation can be expected to be valid when 

Ttf ^ rl/iNan) = {[2rl]/[3N^at,a^])'^' « 1. 

For a trap with wo = Stt x 20 Hz and ro = 20aosc and a 
cloud of iV== 10^atomsof*'^Rbinthe|F = 1,top = -1) 
or \F — 2, mp — 2) state (for which a « 5.45 nm [ll|), we 
have Tts = ri/ro « 0.08 and Ftp ~ 0.09. In excellent 
numerical agreement, using the imaginary time technique 
of Chiofalo et al. [li| and the parameters given above, we 
find the ratio of the kinetic energy to the total energy in 
the ground state to be K/Etot — 0.093. 

Towards obtaining the lowest-energy collective modes 
of a thin spherical shell of superfluid, we consider a trial 
wavefunction for the shell condensate of the form: 
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where A is a dimensionless normalization constant (in 
the thin shell limit) and JF is a smooth real function that 
is negligibly small for \r — Ro\ :s> i?i, (for instance a 
Gaussian: J^{x) = e~^ /^). We note that except for the 
discontinuity at the condensate boundary, the Thomas- 
Fermi wavefunction is in the class of wavefunctions de- 
scribed by Eq.|(21). The function 0(r) has the usual re- 
lation to the velocity of the condensate: v — hV4)/m. 
Rq and Ri are variational parameters corresponding to 
the average radius and characteristic width of the cloud, 
respectively. A description of the cloud in terms of the 
parameters Rq and i?i is expected to capture the salient 
features of the isotropic, £ = 0, collective modes of the 
superfluid shell. 

Collective modes. Starting with the variational wave- 
function in Eq.Q and an appropriate velocity field, v(r), 
standard methods yield equations of motion for the col- 
lective coordinates Rq and Ri (see, e.g., |j|). We start 
by noting that the energy of the condensate for a wave- 
function of the form Eq.Q can be written as: 
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where Ueff is equal to the energy of the cloud if the 
phase (f) does not vary in space and acts as an effective 
potential for the parameters Rq and Ri . It can be written 
in a physically transparent form as a sum of contributions 
from zero-point (or confinement) energy, potential energy 
from the trap, and interaction energy: Ueff[RQ,Ri] = 
Ezp + Etr + Eint- For the trial wavefunction in Eq.®, 
in the thin shell limit (i?i <C Rq), we find 
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to lowest nonvanishing order in Ri/Rq, where c^p, Ctr, 
and Cint are independent of Rq and i?i and are deter- 
mined by the form of the function T in Eq. . 

The variational energy, Eq.© can be used to find two, 
low-energy, collective excitations of the superfluid shell: 
one in which the width, i?i , oscillates around its equilib- 
rium value (the "accordion mode" ) and another in which 
the average radius of the cloud, Rq, oscillates around its 
equilibrium value, ro (the "balloon mode"). For simplic- 
ity, we assume that in the accordion mode, the mean 
radius of the shell stays fixed at tq while the width os- 
cillates and that in the balloon mode the width, i?i, re- 
mains fixed while the mean radius oscillates. While it is 
clear that any exact solution will couple changes in the 
width to changes in the mean radius, the "decoupled" 
oscillations we consider here can be expected to illumi- 
nate the correct low-energy physics. Indeed, in the thin 
shell limit, we find that oscillations in i?i do not affect 
the mean radius Rq and that oscillations in Rq affect 
the width, i?i, at a negligible level for small oscillations 
about equilibrium (smaller by a factor [Rq — ro]/ro). 

The Balloon Mode. For this mode, we consider a ve- 
locity field of the form V(, = l3f, where /? is a variational 
parameter, or equivalently, (piy) — (3mr /h. As a lowest- 
order approximation, we hold Ri fixed at its equilibrium 
value, ri, and only allow Rq to vary in time. By con- 
structing a Lagrangian for the parameters Rq and (3, we 
find /3 = -Ro =^ Vf, = Rqt and an equation of motion for 
Rq: 
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In the thin shell limit, Eq.Q yields R*^ = tq and a 
frequency of small oscillations around this equilibrium 
value: 



Wb ~ wo + 0{rl/Rl). 



(5) 



We note that the thin-shell approximation imposes a con- 
straint on the amplitude of oscillations in Rq: i?™™ ^ ri . 
The Accordion Mode. For this mode, we consider a ve- 
locity field for the condensate of the form v^, = (3{r — rQ)f 
(equivalently: (/i(r) = /3m(r — ro)^/(2fi,)) and allow Ri to 
vary in time while holding Rq fixed at its equilibrium 
value, ro. Following the same procedure as for the bal- 
loon mode, we find (3 = Ri/Ri => Vq = f{r — ro)i?i/i?i 
and an equation of motion for Ri : 
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with mljf = A-KmA^N J'^^q^T'^{q)dq. Using the fact 
that dUeff/dRi = at equilibrium, the frequency of 
small oscillations of i?i about its equilibrium value can 
be written: 
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The thin-shell approximation imposes a constraint on the 
amplitude of oscillations in Ri: i?™°^ ^ ?'o- In the limit 
of weak interactions, Eint <C Etr, Eq.^ reduces to oja = 
2a;o- Thus, in the weak interaction limit, the frequency 
of this mode is equivalent to that of the breathing mode 
of a spherical condensate cloud. In the limit of strong 
interactions, Eint ^ Ezp, we have Eint — '^Etr and find 
Wa = -s/Swo- This result should be compared with the 
strong interaction limit of the spherical breathing mode: 

We note that the modes in the thin-shell limit de- 
scribed by Eq. 13) and Eq. have a structure identical to 
that of a one-dimensional condensate, corresponding to 
its 1-d sloshing and breathing modes, respectively. For 
thicker shells, we expect corrections to our results that 
couple the i?i and Rq degrees of freedom. In fact, in the 
limit that Rq —^ we expect the balloon and accordion 
modes to tend to the breathing mode and next radially 
symmetric mode (n = 2) of a filled spherical condensate 
with Ro — 0. 

Expansion. The dynamics of the spherical shell upon 
release of the trapping potential has noteworthy features 
absent in the case of the filled sphere. Upon release, the 
initial confinement of the condensate causes the outer 
edge to expand outwards and the inner edge to collapse 
inwards. As a result, the system can potentially exhibit 
accumulation of mass at the center, and the condensate 
can interfere with itself when diametrically opposite re- 
gions come together. 

The timescale of expansion can be estimated within 
the thin-shell approximation (where the dynamics of the 
width, i?i, do not affect the mean radius, i?Oi which we 
approximate as fixed at rp). Taking the function J- in 
Eq.© to be a Gaussian and evaluating the different en- 
ergy contributions as before, we find that energy conser- 
vation between the instant the trap is switched off and 
later times t gives the relationship 
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where the time argument of Ri{t) on the r.h.s is sup- 
pressed, and Ri{0) is the characteristic width of the con- 
densate shell before expansion. Assuming the initial en- 
ergy is dominated by the interaction energy, we find that, 
on the time scale for which the shell expands enough 
to become dilute but not enough to reach the center, 
R\{t)/R\{Q) K. 2ujlt^, in contrast to a filled sphere for 
which Rl{t)/Rl{0) « (2/3)wgi2 g^ Hence, for typical 
parameters used in this paper, a shell of initial thickness 
5 /im and radius 50 ^m should expand to a thickness of 
Ri{t) — 20 /im on a time scale of around 20 ms, which is 
amenable to experimental detection. 

To quantify the physics of the expanding shell, we per- 
formed a numerical time evolution of the initial conden- 
sate wavefunction (given by the numerically-obtained re- 
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FIG. 1: Expansion of a thin spherical shell, (a) Evolution 
of the density as a function of radial position and time, (b) 
Snapshots of the density profile at times i = 0, 3, 6, 9, 12 (from 
the lowest curve, subsequent curves are shifted up by 100 for 
clarity) . 



suit discussed earlier) after release from the trap, includ- 
ing interactions. This expansion process was obtained 
using the real-time synchronous Visscher method 9] to 
integrate the time-dependent GP equation: i dip/dt — 
SH/Sip*. The results are shown in Fig.Q] 

Two general features of the expansion deserve discus- 
sion. First, mass accumulation in the center can result in 
a density greater (but not necessarily much greater) than 
the initial density. We note that repulsive interactions 
between the atoms prevent the density from becoming 
as large as it would in the non-interacting case. Sec- 
ondly, interference fringes appear after the inner radius 
of the cloud has reached the origin, demonstrating self- 
interference of the condensate. For two Gaussian con- 
densates initially separated by a distance D, the fringe 
size at long times is given by 5r = 27rHt/{mD) |4| in the 
absence of interactions. The free expansion of an ini- 
tially thin Gaussian shell is straightforward to calculate 
and we find that the fringe spacing at long times is iden- 
tical to that of two Gaussian condensates, but with the 
initial separation, D, replaced by the initial diameter of 
the shell, 2ro. This implies Sr = ntojo a^^^/ro, which for 
To = 20aosc gives a fringe spacing at time tu/Q = 10 of 
i5r — 1.6 Oosc- This compares to an average fringe size ob- 
served in our numerics of about 1.2aosc- The difference 
in precise values presumably results from the effects of 



interparticle interactions and the non-Gaussian shape of 
the initial wavefunction. 

The dynamics of the shell upon trap release has dis- 
tinct parallels with the Bose-nova collapse of a conden- 
sate when the interatomic interactions are switched from 
repulsive to attractive. For the shell, the initial implosion 
is caused by the quantum pressure of the condensate forc- 
ing itself to fill the low-density region at the center. Sim- 
ilar to the Bose-nova case, mass buildup near the origin 
is followed by a relaxation and expansion on timescales 
comparable to the trap frequency. An important ques- 
tion, given these parallels, is whether three-body recom- 
bination and subsequent "loss" of atoms is appreciable 
in the case of the released shell (as it is in the case of 
the Bose-nova collapse [if). These losses are described by 
the equation dn/dt = —K^n^ where K^ is the three-body 
loss rate y. Concentrating on the time in our numerics 
with the largest density (at t ss 6.5/ujo), we estimate that 
the density in the central plateau is n ~ SSO/a^gj. over 
a radius of about aosc- Assuming that this density per- 
sists for the entire time between snapshots, St ss S/cJo, 
and taking K3 = 4.0 x lO""^" cm^/s Q, we find an up- 
per bound on the number of particles lost in this region 
during this time slice of AA^ ~ 0.09, making three-body 
recombination negligible for the case considered here. 

The effects of mass accumulation would be enhanced if 
only the inner edge of the trap were removed, suppressing 
the outward expansion of the condensate. This more dra- 
matic case was considered by Zobay and Garraway ;.Q| , in 
which they modeled an initially shell-shaped trap quickly 
switched into a harmonic trap. Similar features of mass 
accumulation and self-interference fringes were found in 
the case of this bubble trap. For a shell with the param- 
eters given in this paper, we can estimate the timescale 
for collapse and the kinetic energy gain in this scenario 
by considering a small cavity of radius R2 at the center of 
a condensate. In the Thomas-Fermi approximation, the 
dynamics of such a cavity can be mapped to the standard 
hydrodynamics problem of a collapsing bubble in a fluid 
governed by the equation |l2l | 
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Here, p = mn is the condensate mass density, p is the 
pressure at radius r and po = n^g/2 is the pressure far 
from the cavity. Integrating Eq.© at the edge of the 
bubble (r = R2) and making the substitution R2{t) — 
R2{Q)x^/^ gives the time for complete collapse in terms 
of the initial radius, i?2(0): t/ « 0.915i?2(0) ^p/Po- The 
kinetic energy gained by the particles upon reaching the 
center is given by Eke = 47r/3poi?2(0). For a cavity of 
radius 40 /J,m and quantum pressure of magnitude po — 
1 X 10~^^ erg/cm , the collapse time is on the order of 
100 ms and the kinetic energy gained per particle is on 
the order of 1 nK. The small cavity treatment suggests 



that the collapse of the inner radius can be accompanied 
by a measurable amount of kinetic energy gain and mass 
accumulation. 

In conclusion, motivated by possible new trapping po- 
tentials for Bose-Einstein condensed atoms, as well as 
inhomogeneous phases of dilute ultracold bosons in opti- 
cal lattices or in mixtures, we have explored the collective 
modes and expansion dynamics of a superfluid confined 
to a spherical shell. The two breathing modes we find 
offer a possible way of distinguishing a filled spherical 
condensate from a one which is hollow or surrounding 
some other phase. The expansion properties of the shell 
after release from the trap are found to have some notable 
similarities with Bose-nova physics, particularly mass ac- 
cumulation and self-interference. 
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